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Abstract 

The 1-loop computation of the processes — > 1% has been performed at an arbitrary 
cm. energy for any pair of MSSM neutralinos. As an application suitable for Dark 
Matter (DM) searches, the neutralino-neutralino annihilation is studied at the limiting 
case of vanishing relative velocity, describing the present DM distribution in the galactic 
halo; and at a relative velocity of about 0.5, determining the neutralino relic density 
contributions. The most useful situation is obviously for i = j = 1, but the case of 
non-identical neutralinos may also be useful in some corners of the parameter space. Our 
results are contained in the FORTRAN code PLATONdmgZ, applying to any set of real 
MSSM parameters. Numerical results are also presented for a sample of 6 MSSM models, 
describing the various possible neutralino properties. A comparison with other existing 
works is also made. 
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1 Introduction 



The nature of Cold Dark Matter (DM), constituting almost 23% of the energy of the 
Universe, is one of the most exciting subjects of physics today PQI2]- Within the minimal 
R-parity conserving Supersymmetric (SUSY) framework, a most obvious candidate for 
such matter is of course the lightest neutralino Xi CP- A striking signature for such DM 
would then be the detection of 7-rays obtained through the annihilation of two neutralinos 
0- 

The spectrum of most of these 7-rays should be continuous 3j. But once the neces- 
sary sensitivity is reached, sharp 7-rays could also be observed, induced by neutralino- 
neutralino annihilation at rest, to 77, 7Z or a photon together with a neutral Higgs 
particle. Observing such 7-rays, with the predicted ratio of intensities, would be a really 
great discovery. 

For identical annihilating neutralinos in the vanishing relative velocity limit, the pro- 
cesses X1X1 ~^ 77 an d X1X1 ~~ > 1% have already been studied in jlj and j^j respectively. 

Subsequently, x^X°j ~~ y 77 nas a l so been studied in jHj, at any relative velocity and 
for any neutralino pair. This study is based on analytically expressing the amplitudes, 
in terms of Passarino-Veltman (PV) functions Such non- vanishing relative velocity 
results may be useful in estimating the neutralino contribution to the DM density [3]. 
In some, admittedly extreme corners of the parameter space, where the two lightest 
neutralinos may be exactly degenerate, the (i 7^ j)-case may also be useful. 

The same analytic amplitudes were then also used to study the reverse process 77 — > 
X^X°j i n a 77 Linear Collider (LC 77 ) [SJ; while the related amplitudes for gg — > XiX.'j an d 
99 ~^ Xi9i were used to determine the corresponding production of a pair of neutralinos 
0, or a single neutralino [JU], at the CERN LHC. 

FORTRAN codes supplying all 1-loop contributions to these processes, for any set 
of real MSSM parameters at the electroweak scale, and any neutralino pair, may be 
obtained from ^T] . These should be useful for checking the consistency of the neutralino 
DM identification, using collider experiments. The importance of such consistency checks 
can hardly be overemphasized. 

The annihilation process X1X1 ~> 77? f° r the lightest neutralino, at any relative ve- 
locity, has also been studied recently by [T2]. In this reference, an automatic numerical 
method is presented, which directly calculates the needed 1-loop amplitudes, starting from 
the Feynman diagrams. As pointed out by [T2], their results for 77- annihilation, perfectly 
agree with those of [I] and 

However, jT2j also gives results for X1X1 ~ > 1% a t v u — 0-5 and vn = 0. Discrepancies 
appear though, when comparing these results with those of which only exist for 
Vn = 0. An independent calculation of this process seems therefore required. As for the 
77-case, it would be advantageous to have 7Z-results for any pair of neutralinos, at any 
relative velocity, which may allow also the study of the reverse process in a future ej 
linear collider. 

Therefore, in this paper we present an analytic calculation of the 1-loop process x° X°j ~* 
7Z, following the same philosophy as in [HI IE]- The results are based on analytically 
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expressing the helicity amplitudes in terms of PV functions, and they are valid for any set 
of real MSSM parameters 1 . They are contained in the FORTRAN code PLATONdmZ, 
which is herewith released in 

As we will discuss below, the results announced here agree with those of ^2] for i=j=l, 
at vanishing relative velocities, while they deviate from those of [3]. In the vn = 0.5 case 
though, a perfect agreement with ^2] exists only for five of the six considered models, 
while for the other one the agreement is only approximate. Our procedure is presented 
in Section 2 below, the results in Section 3, and an Outlook is given in Section 4. 



2 Procedure 

The process studied is 

Xi(k,Xi) + x j (h,X2)^j(pum)+Z(p 2 ,fM 2 ) , (1) 

where the momenta and helicities of the incoming neutralinos and the outgoing vector 
bosons are indicated explicitly. Generally, the incoming neutralinos may be different. The 
corresponding helicity amplitudes, satisfying the standard Jacob- Wick conventions [T5j . 
are denoted as F Xi x „ 2 (0), where 6 is the cm. scattering angle of the outgoing photon 
with respect to the incoming neutralino Xi- 

According to ^Hj, the x°iX°j antisymmetry due to the fermionic nature of the neutrali- 
nos, obliges the helicity amplitudes to obey 

Fl, rM (d) = (-1)" 2 F{l Xi . M (n - 9) , (2) 
while the CP symmetry for real SUSY-braking parameters implies 

/"a,, a, = (-i)^-^ v^Fl^je) , (3) 

where rji = ±1 is the CP eigenvalue of x° j!4j . Relation (j2J) is very important, since it is 
used below to reduce the number of needed independent diagrams. 

In the calculation we use the 'tHooft-Feynman gauge £iy = £z — £7 = 1, together with 
the standard linear gauge fixing. 

The complete set of needed diagrams consists of: the box diagrams of FigHJ the 
bubbles and the initial and final triangles of Fig|2l the t-channel triangle diagrams of 
FigE3 and the 7 — Z self energy contributions 2 in FigEJ In all cases, the full internal lines 
denote fermionic exchanges, the broken lines scalars, and the wavy ones gauge bosons, 
while the arrowed broken lines denote the usual FP ghosts 3 . The external momenta and 
the polarization vectors of the outgoing gauge bosons are indicated in parentheses, in the 
figures. 

1 As usual, the vacuum velocity of light c, is taken as the unit of velocities. 
2 We thank F. Boudjema for drawing our attention to this contribution 
3 The only diagram of this type appears in the last line of in Fig|5J 



3 



Taking advantage of the Majorana nature of the neutralinos, the direction of the 
fermionic line is always chosen from to Xji so that the wave functions of x° are 
respectively described by positive (negative) energy Dirac solutions. 

We next enumerate the main steps of the calculation. 

We call the first 10 boxes in FigQ " direct" . The contribution of the corresponding 
boxes with 7 Z exchanged, is then determined from them, by enforcing (J2J). Therefore, 
to 

jpij (n\ direct Box ( a\ 

we should add 

to take into account the (7 — Z)-crossing contribution; thereby greatly facilitating the 
computation. 

The last 8 boxes in FigJTJ called " twisted" , have been checked to satisfy (J2J) by them- 
selves. This is also true for the bubble and initial and final triangle diagrams of FigJ21 

The third set of needed diagrams consists of the t-channel triangles of FigEl whose 
contribution we denote as 

Fl x „(er^ . (6) 

To this we should add the contribution of the corresponding u-channel triangles, obtained 
by interchanging x® <-> xp (with the arrows always running from x°i to X?) m FigOl By 
enforcing (J2J), these u-channel triangle contribution is given by 

Fi M .^M^ glc = (-^'^Km^J* - ef-^. ( 7 ) 

Finally, the needed 7 — Z self energy contribution is is shown in FigEJ where the left 
diagram describes the contribution due to s-channel exchange of 4 (H°, h°), while the right 
one describes the contribution induced by a t-channel neutralino exchange. Associated to 
this t-channel self-energy contribution, there exist a corresponding u-channel one obtained 
by enforcing (J2J) on it, through 

Km-^J )^ enersy = (-^'^Hl^uj* - °r scli energy - (8) 

To summarize, the complete helicity amplitudes are given by the sum of the contri- 
butions of the diagrams in FigsHHH together with the contributions appearing in (0E1IE1)- 

As already said, the Aa (9) amplitudes are expressed analytically in terms of 
the PV functions. We have made several tests on these results, in order to eliminate, as 
much as possible, the possibility of errors. These, we enumerate below: 

Requiring (J2J) for the contributions of each of the 8 twisted boxes in FigG] and of the 
diagrams of Fig|21 already provides generally stringent tests of their correctness. 

Moreover, for the 10 direct boxes and the first four twisted boxes in Fig^ which are 
similar to the box diagrams contributing to the 77-amplitudes jHj, we have checked that 

4 This contribution automatically satisfies (J2J). 
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the 7Z diagrams smoothly go to the 77 ones, as p| and the Z couplings are replaced 
by the photon ones. 

The calculation of the 10 direct boxes of Fig^ and of the t-channel triangles of FigEl 
for which no symmetry constraint is available 5 , has been checked several times. 

Finally, we have also checked that our results respect the correct helicity conservation 
(HC) properties at high energy and fixed angles [H . Such tests check stringently the 
mass-independent high energy contributions for both transverse-Z and longitudinal- Z 
amplitudes. Any seemingly innocuous misprint, could not only violate HC, but also 
transform the expected logarithmic energy dependence of the dominant amplitudes at 
high energy, to a linear or quadratic rising with energy, thereby supplying a clear signal 
of error. We have had ample experience of this, during our checks. 

In addition to these, we have, of course, assured that all UV divergences 6 , as well as 
any scale dependence, cancel out exactly in the amplitudes, and that J3j) is respected, for 
real MSSM parameters. 

We next turn to the quantities needed for DM studies of the process (0) . These are 
expressed in terms of the helicity amplitudes as 



where describes the relative velocity of the x°x°-pair, implying 

vl + vfA- 



6 5(mf + rrij) — 12mjmj(mf + m|) + ^Im^rr^p 
+ Vij 8{m t + rnj) 4 



(10) 



up to 0(vfj) terms. A numerical search indicates that for < 0.7, the terms up to 0(vfj) 
in (fTTl|) should be adequate. 

The transverse part of E»j is obtained by discarding the \ii = (longitudinal Z) 
contributions in the helicity summation in JOJ). 



3 Results and comparisons 

For understanding DM observations from neutralino-neutralino annihilation to 7Z in e.g. 
the center of our Galaxy or in nearby galaxies like Draco ^H], the quantity (Q should 

5 For the 77-case, checking the 10 direct boxes was helped by requiring the photon-photon Bose 
symmetry gj], which is not available here. 

6 It is amusing to remark that if A = 1/e — 7 + ln(47r) is replaced by zero, as is done by default in |17j . 
and the dimensional regularization scale is chosen as /idim = %, then the 7 — Z self energy contribution 
vanishes at the 1-loop level. In this case, the diagrams in FigQJmay be ignored. 
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be known at ~ 10~ 3 lj. At so small velocities, the relative orbital angular momentum 
of the x°x°-pair must vanish, and the system must be in either an l S or an 3 Si state. In 
such cases, the angular distribution of dT^j/dcosO is flat. 

Angular momentum conservation implies that the 1 S'o state can only contribute to 
transverse helicities for both 7 and Z\ while the 3 S± state can also give non-vanishing 
longitudinal-Z contributions. 

If the two neutralinos happen to be identical, Fermi statistics only allows the 1 5 , -state 
at vanishing relative velocities, so that 7 and Z, are both transverse. At higher velocities 
though, like e.g. Vu = 0.5, longitudinal-Z amplitudes can also arise. 

On the other hand, if i ^ j, longitudinal-Z amplitudes can also contribute, even for 
vanishing relative velocities. 

We also note that for i = j, the angular structure of dJlij/d cos 9, at non-negligible 
relative velocities, is always forward-backward symmetric. But for i 7^ j, this is not true 
any more. Depending on the content of the two neutralinos, dUij/d cos 9 is sometimes 
peaked in the forward region, and others in the backward; compare For t>y < 0.5, 
dTjij/d cos 9 was found to be flat, in all examples we have considered. 

Next, we turn to the specific features of our approach based on the PV functions, 
whose definition is known to be singular close to threshold and at the forward or backward 
angles 7 [17]. Thus, for relative velocities of ~ 10 -3 , or for angles in the forward or 
backward region, extrapolations must be done. In all examples we have considered, these 
were very smooth, with no suggestion of a possible introduction of errors. 

The herewith released code PLATONdmZ calculates dH^jd cos 9 in fb JTT], for real 
MSSM parameters at the electroweak scale, and fixed values of the relative neutralino 
velocity Vij, and 8 cos 9 using [T7j. For 0.1 < Vij < 0.7 and angles away from the forward 
and backward regions, the direct use of the code usually runs without problems. 

The only known exception appears in cases where the sum of the masses of the two 
annihilating neutralinos happens to be close to the Z-pole. Such additional threshold 
singularities are specific for the XiX^ 1% mode, and they have no counterpart in the 
corresponding 77- and gg calculations of j^j. 

The next step is to compare our work, with that of other authors'. The only preexisting 
results in the literature apply to i=j = 1 for v u = and v\\ = 0.5, presented by [T2] . and 
for Vu = presented by [JJJ. To compare with them, we give in Table 1 the results of the 
PLATONdmZ code, together with those of ^2] and pj. We use the same sample of models 
as in 9 [T2]- These models, whose electroweak scale parameters are indicated in Table 1, 
have been named Sugra, nSugra, higgsino-1, higgsino-2, wino-1 and wino-2 by I n 
Sugra and nSugra, the lightest neutralino (LSP) is a bino. In higgsino-1 and higgsino-2, 
the two lightest neutralinos are almost or exactly degenerate higgsinos. Finally, in wino-1 

7 These singularities are solely due to the mathematical definitions used. They do not have anything 
to do with the physical problem and they do not appear in the total amplitude. 

8 To transform it to the usual DM units of cm 3 /sec we should multiply it by 3 x 1CP 29 . 
9 We also use m z = 0.091187TeV" and s 2 w = 0.2319, as in 12 . and m t = Q.llATeV. 
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Table 1: Results for Ejj = %cr(^°^° — > 7Z) at = and = 0.5, summing over all 
7Z polarizations. The transverse Z contributions are also indicated in parentheses, 
for the relevant cases of either % 7^ j or 7^ 0. Previous results for % = j = 1 from [3] 
at vn = 0, and from |T2| at Vxi = and vn = 0.5, are also compared. Inputs are at the 
electroweak scale using the model sample of with tan/3 = 10, and Af = apart from 
A t = — 0.3TeV for m? = 0.8TeV and A t = for m? = 10TeV; all masses in TeV. 





Sugra 


nSugra 


higgsino-1 


higgsino-2 


wino-1 


wino-2 
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0.2 




0.1 






0.5 




20 


0.5 




20.0 




M 2 


0.4 




0.4 






1.0 




40. 


0.2 




4.0 




A* 


1.0 




1.0 






0.2 




4.0 


1.0 




40.0 




M4 


1.0 




1.0 






1.0 




10.0 


1.0 




10.0 




m / 
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0.8 






0.8 




10.0 


0.8 




10.0 
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jZ) in units of 10 27 cm 3 sec 1 ; 
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2.01 • 


10" 




2.60 • 10- 


-6 




0.224 




0.0266 


12.6 




10.1 
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2.03 • 


10" 


-5 


2.61 • 10- 


-6 




0.219 




0.0220 


11.7 




10.1 







1.42 • 


10- 


-5 


1.79 • 10" 


-6 




0.261 




0.0329 


11.7 




10.1 
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2.45 • 


10" 


-5 


5.28 • 10- 


-6 




0.307 




0.0165 


14.2 
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(5.16 • 10- 
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(0.0165) 


(14.2) 
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-5 
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-6 
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0.0166 
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2.63 • 


10" 


-4 


9.64 • 10" 


-S 
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3.28 • 10^ 


2.57- 10- 


-3 


1.04 • 10- 






(1.94- 


10- 


- 4 ) 


(6.02 • 10- 


-5 


) 


(4.12 • 10" 5 ) 




(1.81 • 10~ 8 ) 


(1.07 • 10- 


- 3 ) 


(8.65 • 10- 


- ? ) 


x?x^ 


4.51 • 


10- 


-3 


4.46 • 10- 


-3 




1.59 • 10~ 2 




4.53 • 10- 4 


0.632 




2.90 






(1.06 ■ 


10- 


- 4 ) 


(1.06 • 10- 


-4 


) 


(9.93 • 10~ 3 ) 




(2.59 • 10~ 7 ) 


(6.25 • 10- 


- 3 ) 


(1.24 • 10- 


6 ) 


Uy = 0.5 


xlxl 


2.78 • 


10" 


-4 


1.02 • 10- 


-i 




9.11 • 10" 4 




1.45 • 10~ 8 


2.45 • 10- 


-3 


9.01 • 10- 


-7 




(2.06 • 


10- 


- 4 ) 


(6.15 • 10- 


-5- 


) 


(3.53 • 10- 5 ) 




(3.98 • 10~ n ) 


(1.05 • 10" 


- 3 ) 
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-3 


4.29 • 10- 


-3' 
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10" 
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-5 


) 


(9.72 • 10" 3 ) 




(2.25 • 10" 7 ) 


(4.39 • 10- 


- 3 ) 


(1.31 • 10- 
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and wino-2, the lightest supersymmetric particle (LSP) is a wino, while the NLSP is a 
bino. 

As seen in Table 1, the PLATONdmZ results for i=j=l and V\\ = 0, perfectly agree 
with those of while they deviate from those of i5J. As expected, Z is completely 
transverse, in this case. 

For vu = 0.5 though, longitudinal-Z contributions, are also possible. Because of this, 
in Table 1, we first give the results for the full 7Z-production, while in parentheses, in 
the next line, the completely transverse jZ production is also indicated. As shown in this 
Table, appreciable longitudinal Z contribution at vu = 0.5, only appears for Sugra. 

At vu = 0.5, important discrepancies between our predictions for £ n and those of 
|12j . only appear for nSUGRA, reaching the 40% level. 

We also note that in nSugra and wino-2, En is very sensitive to the relative velocity 
V\\. This can be inferred from the big difference between the vu = and vu = 0.5 
results, for i — j — 1, in these models. For nSugra this is also elucidated in Table 2. Is 
this sensitivity partly responsible for the discrepancy in Table 1, of the present results, 
with respect to those of 12.? And then, why it does not induce any discrepancy in wino-2? 

Table 2: Sensitivity of En = Vn ■ o\\ to vn in nSugra. 



nSugra 


Vu 


En (lO - * ' ! crri 6 sec' L ) 


0.09 


2.65 • 10" 6 


0.2 


3.03 • 10" 6 


0.3 


3.61 • 10~ 6 


0.5 


5.28 • 10~ 6 



In Table 1 we also give results for £12 and £13 for the above 6 models, at v%j = and 
Vij = 0.5. In parentheses, the purely transverse 7Z contributions are also indicated, which 
are generally not identical to the total cross sections. Longitudinal Z production is often 
important in these cases, and sometimes it even dominates the transverse Z contribution, 
at both shown velocities. 

It is also amusing to notice from Table 1 the sensitivities of £12 and £13 on i>y, as it 
changes from 0.0 to 0.5, in the various models. For e.g. £12, strong sensitivity appears in 
the higgsino-2 and wino-2 cases; while for £n a corresponding phenomenon is observed 
for nSUGRA and wino-2 again. 



4 Conclusions and Outlook 

The neutralinos may be the most abundant particles in the Universe, if they really turn 
out to contribute appreciably to its Dark Matter. They are thus, very interesting objects. 



In addition to this, they are very interesting from the particle physics point of view, since 
their Majorana nature allows them to interact in many more ways, than the ordinary 
(neutral) Dirac fermions. Because of this, detailed studies of their properties, both, in 
astrophysical observations and accelerator experiments are welcomed. 

Through the present paper, an extensive analytical study of the 1-loop neutralino 
amplitudes in any unconstrained minimal supersymmetric model (MSSM) with real pa- 
rameters, has been completed, and the related FORTRAN codes have been released [TTj . 

More explicitly, the DM relevant process 

x-x-^iz, (11) 

has been studied analytically here, for any kinematic configuration, while 

XiX°j - 77 , (12) 

has been presented in jU], following the same spirit. The reverse process 77 — ► x°x , 
which is suitable for a LC 77 collider study, has appeared in jHj; while the LHC production 
processes, containing two or one neutralino, have appeared in [§] and [TUj respectively. 

The formalism of all these 1-loop processes is quite common 10 , while the couplings 
are, of course, different in each case. Thus, the 1-loop Feynman diagrams determining 
the amplitudes for neutralino production at LHC, constitute a subset of those entering 
XaX^ 77 production, which in turn comprise a subset of those needed for x°iX°j ~~ > 

The latter, determines also the "reverse" process 

7 + e T -> X- + X°j + e T , (13) 

where an off-shell intermediate 7 or Z, emitted by the incoming e T -line, interacts with 
another incoming photon, producing a pair of neutralinos. Such a process could be studied 
in a future e T 7 Linear Collider, providing further constraints on neutralinos and DM. It 
is straightforward to get the amplitudes for (JT3J), from those of the process (|ll|h studied 
here. We hope to present results for this in the future. 

It would be really thrilling, if we ever unambiguously identify energetic 7-rays coming 
from Space and being associated to Dark Matter annihilation • R would be even more 
so, if, along with the continuous 7-spectrum, we could also detect the sharp monochro- 
matic photons implied by (fT2llTT|) . But even if this turns out to be the case, the neutralino 
DM interpretation will not be sufficiently convincing, unless detailed accelerator studies 
confirm it. The present work contributes towards this. 
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10 To the LHC studies processes receiving tree level contributions also appear. But the expres- 

sions for them are so simple, that no codes are needed. 
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Figure 1: Box Feynman diagrams needed for XiXj — > F un internal lines denote 
fermionic exchanges, broken lines with no arrows describe scalars, and wavy lines gauge 
bosons. The external momenta and the polarization vectors of the outgoing gauge bosons 
are indicated in parentheses. The direction of the fermionic line is always from xl to x®. 
We call the first 10 boxes "direct". The corresponding boxes with 7 «-> Z exchanged, 
are determined by simply enforcing (J2J), to the above "direct" contribution; see text. The 
remaining 8 boxes, called "twisted", automatically satisfy (J2J). 
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Figure 2: The bubble, initial and final triangle diagrams for XiXj ~^ 1%- The exchanges 
are as in Fig^ but some specific exchanges of W and goldstone bosons G in the bubbles, 
are indicated explicitly. S° denotes h° or H°, except in the middle diagram of the last line 
where it also describes the exchanges of ^4° and G°. S denotes a charged scalar particle 
exchange. The left triangular graph in the last line describes the FP ghost contribution. 
These contributions satisfy (0). 
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Figure 3: t-channel triangles. The exchanges are as in FigQ The u-channel triangle 
contribution is obtained from it through (J7J); see text. 
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Figure 4: The needed 7 — Z self energy diagrams. The left diagram describes the s-channel 
exchange of S° = h°, H°. The right diagram describes the t-channel neutralino exchange 
contribution, while the corresponding u-channel contribution is obtained from it through 
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